Supernova Constraints on Models of Neutrino Dark Energy 
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In this paper we use the recently released Type la Supernova (SNIa) data to constrain the inter- 
actions between the neutrinos and the dark energy scalar fields. In the analysis we take the dark 
energy scalars to be either Quintessence-like or Phantom-like. Our results show the data mildly 
favor a model where the neutrinos couple to a phantom-like dark energy scalar, which implies the 
equation of state of the coupled system behaves like Quintom scenario in the sense of parameter de- 
generacy. We find future observations like SNAP are potentially promising to measure the couplings 
between neutrino and dark energy. 



PACS number(s): 98.80.Es 



I. INTRODUCTION 

Astronomical observations of the Type la Supernova 
(SNIa), Cosmic Microwave Background Radiation(CMB) 
and the Large Scale Structure(LSS) strongly support for 
a concordance model of cosmology where the universe 
is flat and composed of around seventy percent of dark 
energy(DE)P,@,Q. The simplest candidate for dark en- 
ergy seems to be a remnant small cosmological constant. 
However, many physicists are attracted by the idea that 
dark energy is due to a dynamical component, such as a 
canonical scalar field Q, named Quintessence]^. The re- 
cent fits to the SNIa data and CMB etc in the literature 
find that the behavior of dark energy is to great extent 
in consistency with a cosmological constant, however the 
dynamical dark energy scenarios are generally not ruled 
out and in fact one class of models with an equation 
of state(EOS) transiting from below -1 to above -1 as 
the redshift increases, Qumtom0,@ is mildly favored0- 
Being a dynamical component, the scalar field of dark 
energy is expected to interact with the ordinary matters. 
If these interactions exist, it will open up the possibility 
of detecting the dark energy non-gravitationally. 

Recently there have been a lot of interests in the lit- 
erature in studying the possible connections between the 
neutrinos and the dark energ y. M , IE JTfl ITU H^ . IhH 
EE EE E3, EE E3, HH I2E O There seem to be 
at least two reasons which motivate these studies: 1) the 
dark energy scale ~ 1CP 3 ev is smaller than the energy 
scales in particle physics, but interestingly is comparable 
to the neutrino masses; 2) in Quintessence-like models of 
dark energy uiq ~ 1CP 33 eV, which surprisingly is also 
connected to the neutrino masses via a see-saw formula 
rnf,/m p i with m p i the planck mass. 



Is there really any connections between the neutrinos 
and dark energy? Given the arguments above it is quite 
interesting to make such a speculation on this connection. 
If yes, however in terms of the language of the particle 
physics it requires the existence of new dynamics and new 
interactions between the neutrinos and the dark energy 



sector. 

In general for the models of neutrino dark energy or 
interacting dark energy, the lagrangian can be written as 



C 



SM 



Hi 



(1) 



where £„ M is the lagrangian of the standard model (SM) 
describing the physics of the left-handed neutrinos, is 
for the dynamical scalar such as Quintessence. Ci nt in 
(1) is the sector which mediates the interaction between 
the dark energy scalar and the neutrinos. 

At energy much below the electroweak scale, the rele- 
vant lagrangian for the neutrino dark energy is given by 



C v + + M v {4>)vv 



(2) 



where £„ is the kinetic term of the neutrinos. For the 
dark energy scalar part of Cj,, two types of models have 
often been considered with one being the quintessence- 
like and another phantom- like[24) . Thus we introduce a 
factor F in the front of the kinetic term of the scalar 



(3) 



When F = 1 it corresponds to quintessence-like and F — 
— 1 for phantom-like. V(<f>) is the potential of the scalar 
field. 

The last term of Eq.(2) is the scalar field dependent 
mass of the neutrinos which characterizes the interaction 
between the neutrinos and the dark energy scalar. In the 
standard model of particle physics, the neutrino masses 
can be described by a dimension-5 operator 



Lf, = -l L l L HH + h.c, 



(4) 



where / is a scale of new physics beyond the Standard 
Model which generates the B — L violations, Il,H are 
the left-handed lepton and Higgs doublets respectively. 
When the Higgs field gets a vacuum expectation value 
< H >~ v, the left-handed neutrino receives a majorana 
mass m v *~ In Ref.|ll| we considered an interaction 
between the neutrinos and the Quintessence </> 

/3-^-h L l L HH + h.c, (5) 
M pi J 
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where f3 is the coefficient which characterizes the strength 
of the Quintessence interacting with the neutrinos. In 
this scenario the neutrino masses vary during the evolu- 
tion of the universe and we have shown that the neutrino 
mass limits imposed by the baryogenesis are modified. 

The dim-5 operator above is not renormalizablc, which 
in principle can be generated by integrating out the heavy 
particles. For example, in the model of the minimal see- 
saw mechanism for the neutrino masses, 

L = hijiLiNRjH + ^M l3 N c Rl N R3 + h.c. (6) 

where is the mass matrix of the right-handed neutri- 
nos and the Dirac masses of the neutrinos are given by 
mo = hij < H >. Integrating out the heavy right- 
handed neutrinos one will generate a dim-5 operator, 
however as pointed out in Ref.^lJ to have the light neu- 
trino masses varied there are various possibilities, such 
as by coupling the Quintessence field to either the Dirac 
masses or the majorana masses of the right-handed neu- 
trinos or both. 

In Ref.0] we have specifically proposed a model of 
mass varying right-handed neutrinos. In this model the 
right-handed neutrino masses Mi are assumed to be a 

function of the Quintessence scalar Mi(<fi) = M^ M v l . 
Integrating out the right-handed neutrino will generate a 
dimension-5 operator, but for this case the light neutrino 
masses will vary in the following way 



-a~£- 2 
e ' M ?ijl L l L HH 



h.c. 



(7) 



With mass varying right-handed neutrinos given above 
we have in [l3j studied in detail its implication in thermal 
leptogenesis. In Ref.[2l| we have studied the possibility 
of detecting the time variation of the neutrino masses 
with Short Gamma Ray Burst. In Ref. [2^ we have 
discussed the implications of the mass varying neutri- 
nos in the cosmological evolution of the Universe. And 
in Ref. |2^| we argued that neutrinos coupled to Phan- 
tom scalar can provide a scenario of dark energy with 
the equation of state crossing the cosmological constant 
boundary of -1. In this paper we will use the recently 
released SNIa data to constrain the couplings between 
the neutrinos and the dark energy scalar. We will show 
that the current data mildly favor the model where the 
Phantom-like scalar couples to the neutrinos. This pa- 
per is organized as follows: in section II we will present 
the formulation and the results on the constraint on the 
equation of state of the coupled system; In section III we 
study the constraints on the couplings of the neutrinos 
to the scalar fields; Section IV is our conclusion. 



II. FORMULATION AND THE CONSTRAINT 
ON THE EQUATION OF STATE 

In this section we will study in detail the scenario of 
the coupled system of neutrino and scalar field described 



by the lagrangian in (2). The equation of motion of the 
scalar field cj> is given by 



where 



dVi _ dM v IM V 



(8) 



(9) 



is the source term by the interaction between the neutri- 
nos and dark energy, with n v and E being the number 
density and energy of the neutrinos respectively and () 
indicating the thermal average. For relativistic neutri- 
nos, the term is greatly suppressed and the neutri- 
nos and dark energy decouple. For non-relativistic neu- 
trinos, the effective potential of the system is given by 
V ef f(<j>) = V(<f>) + n v M v {4>). In the following, for the 
simplicity we take the 0-depending neutrino mass M u {<j)) 
as M v . 

For this coupled system, the energy for each compo- 
nent does not conserve while the total number of neutri- 
nos is constant. It is easy to get that the energy density 
of neutrinos follows 



p v + 3Hp v = n v M v . 



(10) 



The conservation of energy-momentum tensor of the 
whole system gives that the fluid equation of the scalar 
field 6 is 



+ 3Hp^(l + uty) 



-n„M„ 



fill 



with the corresponding equation of state = P^/p^ = 
2 F ^ V \^\ which stands for the uncoupled equation of 

±F<t> 2 +V(<j>) 

state and F has the same convention as that in Eq. ©. 

From Eqs. I|1U|) and Ijllfl one can easily obtain the 
equation of state of the whole coupled system 



Sir! 



a, + 0, 



(12) 



In the following, we will use the recently released SNIa 
data[l| to constrain this class of models. For the low red- 
shift range covered by the geometric constraints of SNIa 
the equation of state of the coupled system can be simply 
parameterized as u> = Wq + W\Z. In the left panel of Fig. 
1 we plot the parameter space of (wq, w\) constrained by 
the SNIa data. One can see the best fitting values are 
around wq = —1.75 and w\ — 1.9, which gives rise to 
an evolving EOS and transiting from below -1 to above 
-1 as the redshift increases as shown in the right panel 
of Fig.l. In our case for a constant equation of state 
the best fitting value corresponds to w = —1.45. On the 
current SNIa fittings we have assumed the matter en- 
ergy density fraction to be in the range 0.2 < fi m < 0.4, 
which is somewhat optimistic. This will inevitably lead 
to some bias but does not affect the physical picture of 
this paper. On the degeneracies between f2 m and dark 
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energy EOS see e.g. 

H EE El Hi Generically one 
needs to add some priors or combined data analysis to 
extract the behaviors of dynamical dark energy yj, 0|. A 
thorough combined observational constraints on the cou- 
plings needs the modification of the Boltzman code and 
investigations on the full set of parameter space, which is 
extremely time-consuming in computation with the mas- 
sive neutrinos (for a preliminary model-dependent study 
see [1^). In this paper we assume a flat prior of the 
Universe, i.e. f2^ = 0. 
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FIG. 1: The left panel: 2 a SNIa limit on the Dark Energy 
model whose equation of state is id = wo + wiz with prior 
0.2 < fl m < 0.4; The right panel: the evolution of w as 
function of scale factor a with the center value given by the 
left panel. 



From equation 1|12|) . one can see that to have w tran- 
sit from below —1 to above —1 as the redshift increases 
we need < — 1 which implies the scalar field to be a 
Phantom-like. In the next section we will use the SNIa 
data to constrain the coupling of the neutrinos to the 
scalar field. We will point out that the interaction be- 
tween the neutrinos and the phantom field is crucial to 
the Quintom scenario of dark energy. 



III. CONSTRAINT ON THE COUPLING OF 
NEUTRINOS TO THE SCALAR FIELD 

Current experiments have put somewhat stringent con- 
straints on the conventional neutrino mass. Atmosphere 
neutrino oscillations show that there is at least one 
species of neutrino with mass > 0.05 c\ 3Jj. Measure- 
ments of LSS matter power spectrum are very sensitive to 
the total mass of neutrino and a combined analysis from 
WMAP@ and SDSSQ gives an upper bound Y>m v < 1.7 
cV at 95% C.L. A recent analysis has optimistically con- 
strained the upper bound to be 0.42 eV. Cosmological 



combined constraints on £m„ are faced with many pa- 
rameter degeneracies and a recent interesting study [32] 
has shown the degeneracy between dark energy EOS and 
neutrino mass, which can loosen the upper bound of neu- 
trino mass. When neutrino mass is varying with time 
the upper bound on current £m„ would be reasonably 
relaxed and Ref.0] has quoted the value to be around 
3 eV. Neutrino density fraction is related to the neutrino 
mass by 



VL„h 2 = 



92.5eU' 



(13) 



where h is Hubble parameter in units of 100 km s _1 
Mpc -1 . For a quantitative study below on the couplings 
between neutrino and dark energy scalar we will take Q, v 
to be around 0.01, 0.03 and 0.06 as specific examples. 

In the presence of the interaction between the dark en- 
ergy scalar and the neutrinos, the neutrino mass will vary 
as the universe expands. In general different couplings 
give rise to different behaviors of the mass variation. In 
the following we will first give a model independent anal- 
ysis and then go to the detail models. 

From Eas. lfTnjl and (fTTJ) we can get 



Pu = Pv^a 



; My 

M v0 ' 



(14) 



where M v $ is the mass of neutrino at present, a being 
the scale factor and 



1 



P0 = i O0oa" 3(1+ ^ ) (l-— / n v a^ 1+w *Um v ). (15) 

P<t>0 JM„o 

From the Eas. l|14|l and l|15|l we eventually get 
Ptot = po[fi m o a_3 + ^"0 a ~ 3 JY- + ^0 O ^(°)]' ( 16 ) 



//0 



with po being the current critical density and 

1 r M » 

X(a)=a-^ 1+W ^(l / n,a 3{1+w ^dm u ), (17) 



(18) 



and the Friedman equation 

H 2 = H 2 E(a) 2 , 
where 

M 

E(a) 2 = [1W~ 3 + n v0 a~ 3 —^ + n^Xia)}. (19) 

In a flat Universe, the observed luminosity distance of 
SNIa can be expressed as 



d L (a) = T(a)c/(aH ), 



(20) 



where T(a) = J a'~ 2 da' /E(a') and c is the speed of light. 
In the considerations above we have assumed that is 
constant for the simplicity of the discussions. 
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Now we move to the detail model. To have a model 
independent analysis we first expand the neutrino mass 
in the powers of the In a with a being the scale factor. 
Setting u — In a wc have: 

M v = M u0 + Ml(0)u+^M'J(0)u 2 + ... 
M'(0) 

= M vQ {l + \p-u + ...). (21) 

Defining P = ^j^- and keeping the leading order in u 
one obtains 

M„=M„ (H-/3u), (22) 

where M v $ denotes the neutrino mass at the present time 
and P is the parameter characterizing the dependence of 
the neutrino mass on u. Substituting Eq. (1221 into Eq. 
(|16[) and following Eq. 1|17[) . and combing the conserva- 
tion equation of the matter including the baryons and 
the cold dark matter 

Pm + 3H Pm = 0. (23) 

we can get 

ptot = p [n m0 a- 3 + n v0 a- 3 (l + Phi a) + Q^Jf (a)](24) 
with 

X (a) = [1 + ^ (1 - a 3w *)]a~ 3( - 1+w *\ (25) 

the Eq. I|19|l now becomes 
E{af = [f7 m0 a- 3 + ^ a _3 (l + /31na) + ^ X(a)]. (26) 

With a uniform prior 0.2 < O m < 0.4 in Fig. 2 we 
delineate 2<r SNIa constraints on the parameter space 
(P,w,j,) of our model, where the SNIa data are taken 
from the Riess Gold sample 1]J. We can see that a non- 
vanishing P is preferred and the data mildly favor the 
model where the neutrinos interact with the phantom- 
like scalar. The parameter space gets better constrained 
with the increase of Normally we expect that the 
constraints of the parameter P should be | p | < 1 for the 
consideration of the small In a expansion, but the current 
SNIa constraint is very weak which lead to the In a ex- 
pansion of M v in Eq. (1221) loses its significance in some 
sense as a large \P\ is not ruled out and even somewhat 
favored. Moreover, the linear parametrization itself is 
not well defined for a large but positive P, because the 
neutrino becomes massless and even be with a negative 
mass on large redshifts. In our Fig. 2 the hatched denotes 
where P > 0. We should point out that in our numerical 
calculations we used no approximations and just set Eq. 
(|22|l as a parametrization. 

Given the limitation mentioned above and the con- 
ventions in the literature for the parametrization of the 
EOS of dark energy one may consider different types of 
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FIG. 2: 2a SNIa limits on the coupled neutrino-DE model 
with neutrino mass evolving as M v = M v q{1 + ftu), where 
u stands for In a. The dots inside show the best fit values. 
The cyan(light dashed), black(solid) and red(dark dashed) 
contours stand for Q. v = 0.01, 0.03 and 0.06 respectively. The 
hatched area denotes where the neutrino mass might be nega- 
tive on large redshifts and physically forbidden. In the fittings 
we've applied the uniform prior 0.2 < Q, m < 0.4. 
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FIG. 3: 2a SNIa limits on the coupled neutrino-DE model 
with neutrino energy density evolving as p v = ^^p^a - ^ 
where the subscript '0' stands for today's value. The dots 
inside show the best fit values. The cyan(light dashed), 
black(solid) and red(dark dashed) contours stand for Q. v = 
0.01,0.03 and 0.06 respectively. In the fittings we've applied 
the uniform prior 0.2 < Q. m < 0.4. 



5 



0.5 



0.0 ■ 



-0.5 



-1.0 ■ 



-1.5 ■ 



-2.0 






— !i v =0.01 




— — £} =0.01 




n -0.03 




a =0.03 




. . . £j =0.06 




\ . • • • Q =0.06 






\ 


\ 


- 






\ 








\ 



















































0.6 



0.8 a 0.6 0.8 1. 



FIG. 4: Effective EOS for the coupled neutrino-DE system 
for the best-fit values of the two models shown in Fig[2]and 
FigEl The cyan(dashed), black(solid) and red(dotted) lines 
stand for f2„ = 0.01, 0.03 and 0.06 respectively. The left panel 
stands for the model with M v — M„o(l + /3u) and the right 
panel for p v = jf^-p^aT^ . 



parametrization of the neutrino mass. For an example 
we consider here another type of parametrization which 
was firstly invoked for the study of dark energy coupled 
with dark matter in solving the coincidence problem|33fl 



Aa^ 



Pv 



(27) 



where A and £ are constant. A is determined by A = t^-, 



n± 

and the variation of the neutrino mass is d In M„/d In a 
— 7FI • Correspondingly we now get 



E(a) = (1 - n mo 



i Li, 



Sir 



There would be interactions between neutrinos and the 
scalar dark energy when £ ^ — Sw^. £ = 3,W0 = —1 
stands for the conventional non-interacting ACDM cos- 
mology and £ = for the "tracking" neutrino which 
would behave the same as dark energy. As shown in 
Ref.J^] it is not applicable that neutrino and DE can 
enter the tracking regime as early as today, since neu- 
trino takes up a very small density fraction and dark en- 
ergy only comes to dominate the universe very recently. 
Hence we will restrict our attentions on the parameter 
space where £ > 0. 

In Fig. 3 we delineate the 2a SNIa constraints on the 
parameter space of (^Wtp). The cyan(light dashed), 
black(solid) and red(dark dashed) contours stand for 
fl u = 0.01,0.03 and 0.06 respectively. We can see that 
current SNIa constraints are very weak and the data fa- 
vor mildly the model where neutrinos interact with the 



phantom-like scalar. It is noteworthy that here the best 
fitting values show £ + Sw^ > 0, which corresponds to an 
increase of neutrino mass with the redshift, and the be- 
havior is similar to our first example. We find for the two 
parametrizations when the prior of f2 m is relaxed the pa- 
rameter space would be constrained less stringently, but 
a phantom-like DE coupling with neutrinos is still mildly 
favored. 

More recently the authors of Ref.|3fJ made the dis- 
tance measurements to 71 high redshift type la super- 
novae discovered during the first year of the 5-year Su- 
pernova Legacy Survey (SNLS). SNLS will hopefully dis- 
cover around 700 type la supernovas, which is an intrigu- 
ing ongoing project. When performing the current SNLS 
constraints on the neutrino-DE coupling, we found that 
the current SNLS data have not yet been as precise as 
the Riess Gold sample pj and the conclusion above has 
not changed. 

As shown by Eci. (|12[) . for SNIa data only the prefer- 
ence for a Quintom-like behavior of dark energy (where 
the equation of state gets across -1) is fully degener- 
ate with the preference of a coupling between neutri- 
nos and phantom-like dark energy(w^ < —1). In fact 
the fittings in Section II can also be rephrased as cur- 
rent SNIa constraints on the linearly parametrized EOS 
of single component of dark energy, as shown in Ref. [j| . 
To give a more intuitive example we delineate in Fig. 4 
the effective EOS of the full coupled system given by 
the best fit values of Fig|3 and Fig|3| The cyan(light 
dashed), black(solid) and red(dark dashed) lines stand 
for n„ — 0.01,0.03 and 0.06 respectively. The left panel 
stands for the model with M v — M v q(1 + (3u) and the 
right panel for p v — We can see they all 

show a Quintom-like behavior. Actually from Eas. l|26(l 
and 1)28(1 we can easily work out the analytic forms of the 
effective EOS: 



-1 



^{^ a- 3 (-3 + /3-3/31na) 



n*,[-3(B+l)(l + «v)c 
35a- 3 ]}, 



(29) 



where 



and 



D = n„ a- 3 {l + /3liia) + n c i >0 X(a), (30) 



D 



for the coupling form of Eq. I(27|l and 



(31) 



(32) 



for the form of Eq. l|28|l . which could easily get across 
the cosmological constant boundary for the given best fit 
values shown in Fig|2]and Fig|3] 
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The projected satellite SNAP (Supernova / Accelera- 
tion Probe) would be a space based telescope with a one 
square degree field of view with 1 billion pixels. It aims 
to increase the discovery rate for SNIa to about 2,000 per 
year [3^. In the following we will use the simulated three- 
year SNAP data to investigate to what an extent future 
SNAP would be able to detect or rule out the currently 
mildly favored couplings between dark energy and neu- 
trino. In addition we assume a gaussian prior a m = 0.01, 
it is close to future Planck constraints 1 . 

The underlying fiducial model used here is the un- 
coupled ACDM model with = —1,0^ = 0.7 and 
h = 0.7. It corresponds to (3 = in the coupling form of 
Ea.p^l and £ = 3 in the form of Eq.J23- The simulated 
SNIa data distribution for each year is taken from Refs. 
[35L f36j . As we consider 3-year data, the number for each 
bin will be improved by around 3 times. As for the error, 
we follow the ref. (3j| which takes the magnitude disper- 
sion 0.15 and the systematic error <j sys (z) = 0.02 * z/1.7, 
and the whole error for each data is 



<Tmag(Zi) = \ ^sys( z i) + ~ , (33) 

y TT>i 

where ni is the number of supernova in the i'th redshift 
bin. 

In Figure 5 we show the 3-year SNAP constraints on 
the coupled neutrino-DE model given by Ea. (|22|l . The 
red(light dashed), blue(dark dashed) and green(solid) 2a 
contours stand for f2„ = 0.01,0.03 and 0.06 respectively. 
We find the coupling parameter of j3 is restricted to be 
of order unity and the parametrization of In a has made 
its sense when the neutrino density fraction can be order 
of 0.01. 

Correspondingly in Figure 6 we delineate the 3-year 
SNAP constraints on the coupled neutrino-DE model 
given by Eq. (|2"5j) . The red(dashed), blue(dash dotted) 
and green(solid) 2a contours stand for = 0.01,0.03 
and 0.06 respectively. The behavior of this parametriza- 
tion is not symmetric around the fiducial ACDM model, 
and the constraints on the parameters are similar to that 
of the linear parametrization. Thus the future projects 
like SNAP can improve the precision efficiently and is 
potentially able to detect the couplings between neutrino 
and dark energy. 



IV. DISCUSSIONS AND CONCLUSION 

In this paper we have made the first study on the 
constraints on the couplings between neutrinos and the 



1 We found through our simulations the constraints on neutrino- 
DE coupling would be weaker without the prior on matter, where 
the necessity of the prior can be understood from the following 
Figures 5 and 6. 




FIG. 5: Simulated 3-year SNAP constraints on the coupled 
neutrino-DE model with neutrino mass evolving as M„(m) = 
M„(0)(1 + f3u), where u stands for lna. The red(dashed), 
blue(dash dotted) and green(solid) 2a contours stand for 
Q. v — 0.01, 0.03 and 0.06 respectively. The hatched area de- 
notes where the neutrino mass might be negative on large 
redshifts and physically forbidden. The underlying fiducial 
model has been chosen as the uncoupled ACDM model with 
c = 0, w$ = — 1, Qtj, = 0.7 and h — 0.7. We take a gaussian 
prior a m = 0.01. 
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FIG. 6: Simulated 3-year SNAP constraints on the coupled 
neutrino-DE model with neutrino energy density evolving as 
pv = Pq>a~^ , where the subscript '0' stands for today's 
value. The red(dashed), blue(dash dotted) and green(solid) 
2a contours stand for Q. v = 0.1%, 0.01, 0.03 and 0.06 respec- 
tively. The underlying fiducial model has been chosen as the 
uncoupled ACDM model with £ = 3, W4, = — 1, £l<f, = 0.7 and 
h — 0.7. We take a gaussian prior a m — 0.01. 



scalar field of dark energy from current and future obser- 
vations of SNIa. We have shown that the coupled sys- 
tem of neutrino dark energy can resemble the behavior of 
Quintom and is fully degenerate with Quintom in light 
of the geometric observations of SNIa only. A system 
with the phantom-like dark energy coupled to massive 
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neutrino is mildly favored by current SNIa data where 
the neutrino mass is increasing with the redshift. Such 
couplings are promising to be detected by future obser- 
vations. 
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